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ABSTRACT 


The purpose of the present report is to develop certain matrix 
perturbation techniques which can be used in the dynamical analysis 
of structures where the range of numerical values in the matrices is 
extreme or where the nature of the damping matrix requires that complex 
valued eigenvalues and eigenvectors be used. The techniques can be 
advantageously used in a variety of fields such as earthquake engineering, 
ocean engineering, aerospace engineering and other fields concerned 
with the dynamical analysis of large complex structures or systems of 
second order differential equations. 

A number of simple examples are Included to Illustrate the 
techniques. 


vi 
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I. INTRODUCTION 


In the dynamic analysis of large complex structures the numerical 
computations are frequently complicated by the fact that the range of numer- 
ical values of the matrices is much larger than today's computers can 
comfortably handle, In other problems the nature of the damping matrix 
may be such that the structure does not possess classical normal modes and 
in such a case the additional complications of having to deal with complex 
valued eigenvalues and eigenvectors may overtax the capacity of today's 
digital computers. In still other problems the engineer may be interested 
in the effects of small changes in parameters on the response of a large 
complex structure whose response is known for one set of parameters. In 
all these cases it may be possible to employ matrix perturbation techniques 
to good advantage. 

In a previous report (TM 33-484) the author showed how certain 
matrix perturbation techniques could be profitably applied to the problem of 
the design of subsystems in large complex structures. The purpose of the 
present report is to develop certain matrix perturbation techniques which 
can be used to advantage in the class of problems discussed above. The 
techniques developed can be used to advantage in a variety of fields such as 
earthquake engineering, ocean engineering, aerospace engineering and 
other fields which are concerned with the dynamical analysis of large complex 
structures or systems of second-order differential equations. 
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II. PRELIMINARIES 


In this report we shall be concerned with the dynamical analysis of 
large discrete structures, All linear problems in discrete structural 
dynamics may be expressed in the following form 

Mx+Dx + Kx s f’(t)'| 


x(0) = a 


£<0) = b 


d.O) 


If the problem is formulated in an inertial frame, than the NXN matrix M is 
symmetric and positive definite, while if the system is passive the matrices 
D and K are at least symmetric and non-negative definite. The N vector, x, 
may contain both displacements and rotations, while the N vector, f, may 
contain both forces and moments 


A, Reduction to Canonical Form 

Since M, D and K are symmetric and M is positive definite, the 
transformation y=M x, reduces (1, 0) to the canonical form 


i£+ b£+ cx = a(t) 

where 


( 1 . 1 ) 


i _ i <Ti 

B = M 2 DM 2 =B x 

.4 T 

C = M 2 KM = C 
2 L = M"^f_(t) 




X(0) = M 



y(0) = M®b 


( 1 . 2 ) 


B. Classical Normal Modes 

The system (1. 0) is said to possess classical normal modes if and 
only if (iff) it can be reduced to a set of N uncoupled second order systems. 
A necessary and sufficient condition for the existence of classical normal 
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~ .. . j£UUfcMulft*idL . 


- ^ 


modes is that the matrices B and C commute (Caughey, I960; Caughey and 
O’Kelly, 1964), Thus 


BCsCB or M” 1 DM” i K a M’^KM^D 


(1.3) 


In this case there exists (3) an orthogonal matrix T, such that (3) 

T^BT ~ 0 — a diagonal matrix — = 6. . 8. | 

& ij ij i 1 


T^CT = A — a diagonal matrix — A. . s 6, . K \ 
b l] i] i / 


and 


rp r , 1 ,2 ,N"i 
T a l_t_ ,,t_ * > >t_ J 


(1.4) 


where is 1,2. • «N are the eigenvectors of C. Equation (1. 0) reduces to 


12 + ©2 + Azsr(t) \ 

T 

£=T x 

rp V (15) 

z(0) = T T x( 0) , i(0) = T T i(0) f 
r (t) = T T a(t) 

J 

C. Application of Classical Normal Modes 

In many practical problems, the damping matrix D, or its canonical 
equivalent B, is unknown, From extensive testing, however, it is known 
that the modal damping is in the range of one to two percent of critical. In 
this case one frequently assumes that the damping is classical and writes 

9 i= 2X ft© (‘-‘i 

T 

Mow B = T®T , using (1. 6) and the Cayley-Hamilton theorem we may write 

B = C(C) (1.7) 

That is, the matrix B is a matrix function of matrix C. Or expressed in 
terms of M and K 
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If it is assumed that (1,8) is a basic property of the materials used in the 
structure, then (1, 6), (1, 7) and (1, 8) remain unchanged in form, as the mass 
and stiffness of the elements of the structure are changed. This assumption 
results in considerable simplification in the problem. Using (1, 6) and 
recalling that © and A are diagonal, (1. 5) may be written 


». + 2 X.z, ar.(t) ^ 
X 1 1 XX x' ' 


z i (0)=z i ' h {0)=i i 


i=i, 2* . . N 




Thus 


(1.9) 


£ 

z.(t) S U.(t)z? + v.(t)z?+ f V,(t-T)r.(r)dt i=!,2...N (1,10) 

1 1 1 X 1 u j> 1 


where 


u. (t) s e 
x v 


\h.t r 

x 


i \Fc. i 

cos*XFt+ — -- x - 1 ' sinX*t 

1 If 1 


v i (t) = 


XFC-t 

1 1 i 

— x — sinT^t 


\=VK> 


/ 


(i. n> 


Equations (1. 10) can be expressed compactly in matrix notation as follows 

,t 
“0 


z_{t) = u(t)z.(0) + v(t)zJO) + J v(t-T) £(t) dT 


where 




u(t) ij =6..u. ( t) ; v(t) ij =6 i .v i (t) 


(1. 12) 
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Using Equations {] . 4} and (1, 5) 


^(t)^Tz(t) = Tu(t)T T x(0) + Tv(t)T T i(0)+ f Tv(t-T)T T q(T) dT (1. 13) 

~ 0 


t 

x(t) =Qu(t)Q T Ma + Qv(t)Q T Mb + J Qv(t-T)Q T £.(T) dT 


Using Equation (I. 2), the solution to Equations (1. 0) can be expressed as 

*fc 
0 

where Q is the congruence transformation 

Q = M“^T 

which has the properties that 


(1. 14) 


(1. 15) 


1) q t mq = i \ 


2) Q DQ = © } 
.T 


( 1 . 16 ] 


3) Q KQ = A 

t 

D, Nonclassical Normal Modes 

If in Equation (1. 1) the matrices B and C do not commute, then in 
general it is impossible to reduce (1, 1) to a set of uncoupled second order 
equations. In this case we rewrite Equation (l t 1) in the form 

dw 


w{0) = 




where 


w(t) = 


M^b 


i(t) 


(1. 17) 


As 


■o I ' 
-C -B 


(1. 18) 
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a, 18) 

cont. 


p(t) 


r 


The solution of Equation (1, 17) is 

w(t) = o At w(0')+ f C A(t - T) J ,(T)d T a. 19) 

- ->0 

Alternatively, there always exists a non- singular matrix S, such that 

S ” 1 AS s J (1,20) 

If the matrix A has a full complement of ordinary eigenvectors, then J is 
■a diagonal matrix whose elements are the eigenvalues of A. Furthermore, 
the matrix S has, as its columns, the eigenvectors of A, If A does not 
possess a full complement of ordinary eigenvectors, then J is a Jordan 
matrix whose diagonal elements are the eigenvalues of A, In this case the 
matrix S has, as its columns, the ordinary and generalized eigenvectors 
of A, 

I?j ^he case where A has a full complement of ordinary eigenvectors 

Equation (1, 19) can be written in a more convenient form for computation, 

-1 - 1 
Since S AS = J = A a diagonal matrix A-SAS 


rx lt 

e 


At 

e s s 


V 


X-t 
e j 


-1 


( 1 . 21 ) 


w(tj=S 


f\ V 

_ 1 , 

-it 

\ (t-T) 

e 

S w O H- 

S 

e 1 

l \ - 

~ *1 

D 

L \ J 


S _ 1 £(T)dT 


( 1 . 22 ) 
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Ill, MATRIX PERTURBATION THEORY 


A , Linear Eigenvalue Problems 

Case 1 [’the eigenvalues of distinct! 

Consider the eigenvalue problem 

U.y 1 * [C° + €C l +6 2 C 2 -.1i l aO (2, 0) 


where G*, C^, etc, , are symmetric NXN matrices and G^ is positive 
definite with distinct roots, and € is a small parameter, Let 

( 2 , 1 ) 

where T is the orthogonal matrix which diagonalizes Thus 

Cu.I + A 0 +6C 1 -i-6 2 C*t. - 0^ = 0 (2.2) 


where s T^C^T s Qp^ t Aq a diagonal matrix with elements X^, Since 

P 

C=rA 0 + Y 
j-1 

is a symmetric matrix, we know that there exists an orthogonal matrix Z 
which will diagonalize C, In addition, we know that the columns of Z are 
the eigenvectors cp 1 of C, and both they and the eigenvalues u. , i=l,2. . . , 
are analytic functions of 6. Knowing this we expand u, and cp x in a power 
series in €, Thus 


flNs&+%ei + € 2 s& 


14 . =m n + + 6 m, 


( 2 . 3 ) 


x ‘ iO ’ ^il ' " *i?, 

Substituting (2,3) into (2,2) and equating coefficients of like powers 
of € to zero, we have 
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C Cu i0 I+ A 0 l2o c0 

€ ^iO I+A 0^2 = ” cl ~r C ^0 wU il^r u i2^0 


(2,4) 


(2, 5) 


(2.6) 


From (2, 4) 


^iO c " ^iO i=;:L ' 2 '* ,N 


j’gise. — a column vector whose elements are 
1 all zero except the element 

GMMj 

(^.A^)=-X 10 5.. 

1 

The rg^'s form a complete orthonormal set. From (2, 5) 

Cn i0 I+ A o^i = ££o* u il^O 


(2,7) 


(2.8) 


Since the eg^’s form a complete orthornormal set they span the N space, 
thus any vector in that N space may be represented as a linear combination 
of the cp* ' s. Let 


4 = l a ii^ 


(2.9) 


Substituting (2, 9) into (2. 8) we have 


j=l 


( 2 , 10 ) 


If the inner product of (2. 10) is taken with cg^ and use is made of (2. 7), then 


8 
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“ki ^ \o " ^kO ^ = ” (^o ' C ^ ) ‘ u i 1 6 ik 


(2. 11) 


Since \ Q *\ 0 , i*K then 


n 1 

Ik 


\i = j ,. 


iO kO' i^k 

We note in passing that since C* is a symmetric matrix cAs;C^. 


( 2 . 12 ) 


Hence 


C M 


If i=k, then we have 


a ik \o“ X iO ^ 


( \ _1 i\ n 1 


(2, 13) 


(2. 14) 


To determine ol .. we make use of the normalization properties of the cjf's 


(eg 1 , jg 1 ) = 1 =(4>-4>) + 2e(t£ l 0 , 4 )+ O (€ 2 : 


Equating coefficients of like powers of €, we have 


1 


= a. . = 0 

ii 


Thus 


u. = x. 0 -^ i + o(e 2 ) 


N r lj 

1 i /- \ ii~H) , r^ir^\ 

2 = 2o ■ s 2 j rr^rr + ° (€ > 


m. 

“j "iCT X jO 

Using (2, 18) we may define a matrix Z, where 

Z = 1+ €S 


(2. 15) 

( 2 . 16 ) 

(2. 17) 
(2. 18) 


(2. 19) 


where 
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0 

a l2 a 13 ‘ ' 

' * a lN 

a Z 1 

(f 

0 a 23 * ' 

' ’ 

c 'm 

%2 

. . o 


The matrix Z ha 3 the following properties 
i) Z T Z«I+0(€ 2 ) 

ii) Z T [a 0 +€C +€ 2 C 2 +-"]z = A 0 +€A 1 +O(€ 2 ) 


where 





( 2 , 20 ) 


( 2 , 21 ) 


( 2 , 22 ) 


(2,23) 


(2. 24) 

where are the eigenvalues of and i— 1,2, • * *N are given by (2. 14), 

Case 2 [some eigenvalues of repeated] 

Suppose that the first r eigenvalues of are repeated. We know 
that the corresponding eigenvectors are not unique and hence the orthogonal 
transformation which diagonalizes is not unique, Let 


V 


u 


11 


u 


21 




N1 


* 

T = T 


r t 

r 

0' 

0 

T 


(2. 25) 


where T is any orthogonal matrix which diagonalizes C® and T r is any rXr 
orthogonal matrix. Then 


t *t t # 1 


(2.26) 


10 
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and 


j.. nr n 

fjiv J* £ v rp 


"x I 

10 r 

1 


A^‘ r 

0 J 


« A, 


(2.27) 


whore is repeated r times. 

Properties (2.26) and (2.27) give us the necessary freedom to solve 
the perturbations in the case of repeated roots, The first step is to apply 
the transformation 


£"Tz (2,28) 

where T is any orthogonal matrix which diagonalizes C^, Thus 

[u.I + A Q + 6C 1 + e 2 C 2 ]cpV0 (2, 29) 

Direct application of the matrix perturbation technique, applied to Case 1, 
fails unless (3^ = 0, i^j, i, jt(l,r), In general these conditions are not 
satisfied, however we can always construct a new transformation T* such 
that this condition does hold. 

Let B r =[C?.. ] i, j € (l,r), and let T r be the rXr matrix whose columns 

-J 

are the eigenvectors of 


[|il+ B r ]t r = 0 


(2. 30) 


The new orthogonal transformation (2, 25) using this particular matrix T 
has the property that 

1 * 

T ■' C T = C 


where 0^ is a symmetric matrix with the property that 


C.. = C..6.. i, j 6(1, r ) 

ij 11 lj > J \ / 


(2. 31) 


(2. 32) 


Thus if T* is used in place of T, T^C ^*=0^ has the necessary properties 
so that the matrix perturbation technique can be used. 
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If in Equation (2, I) T r,> is used in place of T, Equation {2, 2) becomes 
ru.I + A n + €C 1 +€ 2 C 2 . * *1 fp 1 ■= 0 (2.33) 

u 1 u J— 

where C^T* T C J T* j = l,2 .... As before, let 


I i 4 2 i 

C£ =: CPq + + 6 C£^ + * • • 




(2. 34) 


u. = n. 0 +€u.i + e u. 2 + ... 

✓ 

Substituting (2. 34) into (2, 33) and equating coefficients of like powers of 6 to 
aero, we have 

6° [u io I+A o ]4 = 0 


[^o I+ V^ = = gl 2o^u4 


[V + V4 = -cM- s ^u*r^ 


etc. 

From (2. 35), we have 


^iO “ ” \o 


2o=^d 


i= 1,2 • • • N ''j 




(^■£4) =6 ij 

G^O’ A o2o) = "\o 8 ij 

J 

The cp^ 1 s form a complete orthonormal set. From (2. 36) 

UV+V^ = " C ^0"^i0^0 


(2.35) 
(2, 36) 
(2. 37) 


(2. 38) 


(2. 39) 


As before, let cp^ be represented as a linear combination of the cp^'s, Thus 
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(2.40) 


N 


2l = > a ji2o 

j=l 


k 


Substituting (2. 40) into (2, 39) and taking the inner product with cg^, we have 

\i(\o - x k0 ) = - (s£ - £1 4 ) - »iA 12 • 4 1 1 

Thus i^k and either i or k£(l,r), we have as before 


If irk, then 


„ - "ik 

\o“\to 


u n = -(4^) s “ s u 


(2,42) 


(2.43) 


If i^k, but i and k€ (l,r), we cannot determine co^,. As before, normalization 
gives us CLjjS0,\/i. 

2 

To determine i^k, i and k€(l,r) we proceed to the 0(£ ). From 
(2. 37), we have 

i #1 i cs2 i 


^i0 I+A 0*^-2 = ^0" U il-r^2 CP 0 


(2.44) 


Since the cg^j's span the N space, let eg* be represented by a linear combina- 
tion of the cg^'s. Thus 

N 

(2 - 45) 


j=l 


JK 

Substituting (2. 45) into (2. 44) and taking the inner product with cg_, we have 

- (aS* s Vi) - Ofr - ^ 8 ik (2 - 46) 


Making use of (2. 40), we have 
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(2,47) 


N 


^ X i0"*k0^ki = “Is a ji^jk ” ^ik " U il\i " U i2 6 ik 


Jk 

M 

if issk, but i and k€(l,r), Then, since 3^=0^ 6^ for j and k€(l,r), 


gives us 


Therefore 


N 


0= (u kr U il )0t kiT- a «®iv"C?L 


j=l+r 


*** • w " w» 1 

ji jk lk 


N 


I «A + & 

„ _ jal+r 


a. 


ik 


Hence 


n. = x^-effi+ofe 2, 
'l iO u 


N 

cp^^y ct.^+o^ 2 ) 

j=l 

where the a., are given by (2.42) and (2,49). 

Using (2. 51) we may define a matrix Z, where 

Z=I+ €S 

where 


rn 

S = -S = 


0 a l2 a l3"' a lN 

°2l 0 a 23 




(2,47) 


(2,48) 


(2.49) 


(2, 50) 


(2.51) 


(2. 52) 


(2. 53) 
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The matrix Z has the following properties, as before 


where 


i) Z T Z = I+0(€ 2 ) 

ii) Z T [a 0 +6C 1 +€ 2 C 2 + * • *]z = Aq+SAj+O^ 2 ) 


10 


^0 = 


h\ 


NO 


\o~No i=l»2, ♦ ••r 




11 


Ai = 


M 


21 


' U 


N1 


(2. 54) 
(2. 55) 


(2.56) 


(2. 57) 


where A^, i= 1, 2 . • . N are the eigenvalues of and the are given by (2, 43). 

r o *i 

Case 3 [some eigenvalues of C close to one another] 

Suppose that the first r roots of Aq are close such that 


I V X jo^ k€ (i,j)€(i,*) 


(2.58) 


In this case we observe that €cg* is no longer guaranteed to be of 0(€). This 
difficulty is easily remedied by writing our eigenequation in the form 



A 0 + 


SC 1 + € 2 c 2 


C£ 


= 0 


(2. 59) 


where 



u 

1—1 

r-( 



,N~r 

A o J 




(2.60) 
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C 1 = C 1 + i 


Vi 

t 

~ 


0 

\ 




o 

f 

0_ 

/ 


( 2 . 60 ) 

cont, 


Thus the problem has been reduced to one in vyhich the first r eigenvalues of 
C° are repeated. In this case the methods used in Case 2 are applicable, 


B, ’’Nonlinear 11 Eigenvalue Problems 

In the case of nonclassicai normal modes we have to deal with an 
eigenvalue problem in which the eigenvalue occurs quadratically, viz. 


[y^I+€Y.(B 1 +€B 2 +, ..J+cjj/ 1 = 0 (2.61) 

* 

where B , i=l,2> • • and C are symmetric matrices which do not commute. 
Let 


6 1 = Ttp 1 


(2.62) 


where T is the orthogonal matrix which diagonalizes C, Substituting (2.62) 
in (2. 61) 

[\?I+ Y^^+^ft 2 * • ■)+ Aq]^ = 0 (2.63) 


where 


T T 

A 0 = T 1 C = 


n\ 


w. 


w. 


• 2 
w 

n 


1 T i 

a 1 = T B T i=l,2 


(2.64) 


16 
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Equation (2, 63) may also be written in the form 


where 


Let 


Substituting 
to zero, we 


From (2,68) 


[Y.H-A] 


(2. 65) 


A 


‘ 0 

1 

A o 

€b 4€ 2 ft 2 


and 


C£ 

Lve'J 


( 2 , 66 ) 


i i . /--.i i r*2 i 

eg =f£ 0 +%+c egg - 


Y. = Y. n +€y. n +€ 
i i0 xO i2 


(2. 67) 


(2, 67) into (2, 63) and equating coefficients of like powers of € 
have 

.0 [\,2 


,i+A„]a4 » 


0 


e [ Y icr 1 "o 

gl [ Y io I + A o]4 = -[ 2Y io Y uSo +Y io' 3l 4] 


we have 


Y i0 = J w i 


\ 


Y i0 = -j w i = Y i0 

^ = e. (real) 

($0>4) = 6 ii 




( 2 . 68 ) 

(2. 69) 


(2. 70) 
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A o~o) 


2 . 

w. 5 . . 

1 ijj 

The s form a complete orthonormal set, 


(2,70) 

eont. 


For each the? e are two § ll s, one corresponding to vt^, the other 


to V; Q , Thus 


§ i+ 

1 q 


" 4 

$ 1 " r 

r i n 

So _ 

ns? 

> 1 o 

- Y i0 


(2, 71) 


We see that 


■^0 *0 


(2. 72) 


t I 4 

Since the cp^ s form an orthonormal set, let us express cp^ as a linear 
combination of the cg^'s, j=l,2 * . . N, Thus 

N 

$1 = 1“^ < 2 ’ 73 > 

j=i 

Substituting (2,73) into (2,69) and taking the inner product with we have 

[ Y i0 - foKi = - 2 4o Y il 6 ik + to < (2 ' 74) 

We shall consider only the case where the Y^ 1 s are distinct; the 
case of repeated and close roots can be handled in the same way as in the 
linear eigenvalue problem. Thus if i^k 


* _ _ T i0 ik 
a ki “ y:b2 <^2 
Y i0 " T k0 

. ± ±j w i B ik 

‘ * “ki " 2 2 

w r w k 


(2. 75a) 


(2. 75b) 


If isk then 
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(2,76) 


Y il 2 ii 


The OC^'b can be shown to be zero by using the normalization condition on 
thec£ l 's. Thus 

Y? = ijw.- f&J. + O^ 2 ) (2.77) 


N ip 1 

m i* i . , V • i Ik k , n ,,2 . 

cp =(^±6^] 2 22o + 0 ^1 

, , w. -w. 

k=l i k 

Using (2,78) we can construct the two NxN matrices 


<7 

"1 “ 


and 


Z„ = Z. 
2 1 


Let us define the matrices S, A^ , 


Q_ — 
0 


w. 


w. 


w 


N 


(2.78) 


1 

€a l2 

ea{ 3 . . 

• ’ 6a lN 


6a 21 

1 

Sa 23 

Sc 4n 

(2.79) 

j4ti 

€a N2 


1 



(2.80) 


(2. 81) 


S = 


a 


12 


& 21 
”T 2 

w 2“ w l 

«N1 

“2 Z 

w N -Wi 


2 2 
w r w 2 


0 • 


H 


IN 


w r w N 


2 2 

15 ^ 
,5 2N 

2 2 
w« -w 
2 N 


(2. 82) 
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The matrix S is skew symmetric, I, e, , 


S T B -S 


r»!i 


?®2Z 


1 A 

2 %N 


Then the matrix Z ^ can be expressed in the form 

Zj = [i+jesn 0 ] 

The nonsinguiar 2Nx2N matrix, T, whoso columns are the 
vectors of matrix A is given by 


T = 


where 


‘ Z 1 

JLl 

, z l n 

z 75 


A ~ j° 0 ~ 

Matrix T can be expressed in terms of S, 0 Q and A^. Thus 


T = 


i+j6sn Q 

i-j<=sn 0 

i +sc? o) 


-jn 0 - £ 

( A i +sn o). 


+ Og 2 ) 


It may be shown, after some tedious algebra, that 


1 1 

"i+jefio 1 

A j +S Oq 

-jOg'*',. 6S 

“ 2 

J-jen" 1 

MJ 

jn-'-es 


+ 0 (£ 


The matrix T has the property that 


(2.83) 

(2.84) 

(2.85) 
eigen- 

( 2 . 86 ) 

(2.87) 

( 2 . 88 ) 

(2.89) 
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+ 0(£ 2 ) (2,90) 

Thus T is the required similarity transformation matrix which will diagona- 
lize A to 0(€ 2 ), 



IV, APPLICATIONS OF MATRIX PERTURBATION THEORY 

Two examples will now be given to illustrate the use of matrix 
perturbation theory, 


A. 


The Matins of a Small Complex Structure to a Large 



Consider the problem of mating a small, but complex, structure 
with M degrees oi freedom to a much larger structure having N degrees of 
freedom, A good example of such a problem is the mating of a space craft 
to its booster. Suppose that the eigenvalues and eigenvectors are known for 
each structure separ&sely, Can one compute the eigenvalues and eigenvectors 
for the composite structure using this information? As we shall now show this 
can be done without having to solve the eigenvalue problem for the (N+M) 
degree of freedom system. 
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is the inertia matrix for the composite structure and consists of the two 

2 

submatriees Mj and 6 is an NXN symmetric positive definite 

2 

matrix and represents the inertia matrix for the large structure. 6 M 2 

is an MXM symmetric positive definite matrix and represents the inertia 

2 

matrix of the small structure, It is purposely written in the form 6 M^ 
to emphasize that the inertial elements of the small structures are very 
much smaller than those for the large structure, Whore 


Ks 


K 1 

0” 

+ g 2 

r K zi 

K 22l 

0 

0- 

T C 

- K 23 

K 2 J 


(3,2) 


is the stiffness matrix for the composite structure and consists of the five 

2 

submatrices K ]( € K^. i=l,2,3,4, 

Kj is an NXN symmetric nonnegative definite matrix and represents 

2 

the stiffness matrix of the large structure and € » i=l,2,3,4 are coupling 

2 

and stiffness matrices for the small structure, € is an NXN symmetric 

2 

nonnegativ x definite matrix, € K22 an NXM symmetric nonnegative 

2 

definite matrix, € K ^ is an MXN symmetric nonnegative definite matrix. 

2 

€ K24 is an MXM symmetric nonnegative definite matrix. These matrices 

2 

are written in the form € K^* i=l,2,3,4 to emphasize that the elements of 
the stiffness and coupling matrices for the small structure are very much 
smaller than those for the large structure. It should be noted that the eigen- 
values for the two systems separately can be of the same order of magnitude 
since they are governed by the eigenequations 

[XjMj + Kj 3 _Xj = 0 


.2 


cix 2 m z *k 2 ]x z = o 




and are independent of 6, 
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(3. 5) 


In Equation (3. 0) D is assumed to be an (M+N)X(M+N) nonnegative 
matrix which satisfies the conditions for classical normal modes, i, e. 

M* 1 KM" 1 D = M" 1 DM* 'K (3.4) 

In Equation (3. 0) the force vector _f(t) is given by 

P(t)l N 
0 JM 

where p(t) is an N vector and acts only on the large structure. 

Consider the eigenvalue problem associated with (3. 0) 

Cli.M+K]x. = 0 

l —i 

[Note: Since by assumption D is such that (3. 0) has classical normal modes, 
the eigenequation (3. 6) does not include damping.] The transformation 

x.sM"^ 2 ^ 1 (3,7) 

reduces {3. 6) to canonical form. 


(3.6) 


[u i I + C°+€C 1 +€ 2 C 2 ]ji 1 = 0 


(3.8) 


where 


r,,- 1/2__ . .-1/2 

i 

o 

i 

o 

^.-1/2 ,,-1/2 
M 2 K 24 M 2 J 


(3.9) 


c 1 - 


2 

C = 


1 

o 

. .-1/2 . ,-1/2 
M 1 K 22 M 2 

[m- 1/2 k 23 m- 1/2 1 
m- 1/2 k 21 m- 1/z 

0 

0 

( 

o 

0_ 


(3. 10) 


(3. 11) 


- 1/2 - 1/2 

If Tj is the orthogonal transformation which diagonalizes 

- 1/2 - 1/2 

and is the orthogonal transformation which diagonalizes , 


let 
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(3. 12) 



The transformation 


£ = Ttp 1 


(3. 13) 


reduces (3. 8) to 

[m.I+Aq+^+^C 2 ]^ 


(3. 14) 


where 


where 


and 


A 0 ~ 


A 01 “ 


A 01 

0 ' 

0 

A 02. 

\o 



f 


X^q , are the eigenvalues of the large structure 


{3. 15) 


(3. 16) 


,M 

1 0 


02 


* X M 


(3. 17) 


MO 


\f Q , are the eigenvalues of the small structure 


1 T 1 
C = T C T 


(3. 18) 


and 


2 T 2 

C =T C T (3. 19) 

The eigenequation (3. 6) has now been reduced to the standard form 
where the techniques of Section 2 may be applied. The exact details will be 
reflected in two numerical examples. 
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By perturbation analysis we obtain 


Z=I*(=S (3,20) 

and 

A=A 0 +^A 1 (3.21) 

where • > is given by (2. 20) and the elements of A by (2, 18). The orthogonal 
transformation which diagonalizes (3. 8) is given by 


T 2 = T Z = T + S TS + 0(€ 2 ) 


( 3 . 22 ) 


The congruence transformation which diagonaLizes (3, 0) is given by 



(3. 23) 


Q 2~~ 


Mj 1/2 T 1 (I+€S 1 


M, 


T 2 S 3 


eM- 1 / 2 Tl s 2 


k — N modes 




t 2 (i+6S 4 ; 
~ M modes - 




where S^, S£, S^, S^_ are the partitions of S 


(3. 24) 


S = 



N 

M 


(3.25) 


It will be seen from (3. 24) that the first N modes are ordinary ,/r global 
modes. By that we mean that they are simply extensions of the modes of 
the main structure into the secondary structure. The next M modes are 
however "local modes". By this we mean that the response of the secondary 
structure is O(^f), or 0 (^ 2 ) times as i ar g e as that of the main structure. 
That is the motion is highly localized in the secondary structure. 

Consider now the forced vibration problem (3. 0) with a = _b = 0. Using 
Equation (1. 14), the solution of (3,0) is given by 
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x(t) = qJ v(t-T)Q T f_(T)dT 


where v(t') is a diagonal matrix whose elements are given by (1, 11). 
Defining the partitions Ck, 1=1,2, 3 of Q by 

N M 

„ r°i i €q 2 i n 


where Ck, i=l,2,3,4 are 0(1), 


Q 3 I € Q 4j M 


Defining the partitions x^, i=l,2 of x by 


x(t) = 


x^ N 


x 2 /M 


Defining the partitions v^, i=l,2 of v by 


(3.27) 


(3.28) 


v(t) = 


N M 

" V 1 I °] N 


0 v 2 M 


(3.29) 


have 


Introducing the partitions (3.27), (3.28) and (3.29) into (3.26) we 


“ (t) _ 2 1* )j [Q 4 v 2 (t-T)Qj + Q 3 v 1 (t-T)Q^]£(T) 


(3. 30; 


Thus, we see that the response of the main structure and the 
secondary structure depend only on the frequencies in the input and the 
distribution of the eigenvalues. In particular, the response of the secondary 
system is of 0(1) with respect to € and not o(^) as one niight have expected 
from the eigenvectors of the local modes. In order to make the exact details 
of the procedure clear we shall now consider two numerical examples. 
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Example AX (Distinct Eigenvalues) 


In Equation (3. 0) Let the matrices M and K be given by 

1 


Mss 


1 


K - 


“2 

-1 

o' 

H*2€ 2 

“o 

0 

0 

-1 

2 

0 

0 

0 

-1 

0 

0 

OJ 


u0 

-1 

0 


M 


- 1/2 


1 


1 


L *?J 


(3.31) 


(3.32) 


(3. 33) 


.- 1/2 


Using the transformation x = M j; reduces (3. 6) to the canonical form 

(3. 34) 



r 2 -i oi 


"0 0 0' 


m+ 

-12 0 

+ € 

0 0-2 


1 

- 0 0 2. 


.0 2 0. 



(l? 1 = 0 


Now, the matrix 


2 .<1 

• 1 2 


is diagonalized by the orthogonal matrix 

1 1 ” 


✓z 

/2 

1 

1 

/2 

"v/2 


and the matrix [2] is diagonalized by the matrix [l]. Thus 

"l 1 


T = 


Jz /I 

1 1 

/2 "/2 
0 0 1 


The transformation J[_ 1 =Tcp 1 reduces (3. 34) to the standard form 


- 

r 1 i 


r 0 

0 

-v/2” "] 

u.x + 

3 

+ € 

0 

0 

/r 

i 

L 2. 


Lvz- 

n 

0 J 




(3. 35) 


(3. 36) 
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This is a case of distinct eigenvalues (1,2, 3) and hence the standard matrix 
perturbation method can be used. Thus 


z = i+es = 


T 2 = TZ 



1 

0 

-e/2] 

= 

0 

1 

-e/2- 


e/2* 

La# 

£A 

1 

j 


~ J_ 


-2e~ 



A 

A 



_ 

1 

1 

0 



v/2 

"a 



e/2 

e/2 

1 



(3. 37) 


(3. 38) 


Q, 


M 


- 1/2 


•_L 

J- _■ 

A 

A 


__i_ 

A 

~A 

A 

A £ 


0 

-1 


(3. 39) 


has the following properties 

i) QjMQ 2 = l+0(t 2 ) 

ii) Q 2 KQ 2 = A q + €Aj_ + 0(t 2 ) = 


+ 0(c ) 


(3.40) 

(3.41) 


Thus we see that the first body, or global mode, has eigenvalue 1 

and eigenvector (-~,-^“ > /2j and is simply an extension of the first mode of 
S A v /2 

the main structure into the secondary structure. The second body, or 
global mode, has eigenvalue 3 and eigenvector and is simply 

A, A 

an extension of the second mode of the main structure into the secondary 
structure. 

The third mode is a local mode and has the eigenvector [-2€, 0, €~ ]. 
For example, if e = 10~\ then the motion of the secondary structure is 
fifty times as large as that of the main structure. 

If, for example, we take the matrix D to be 
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Then 




-1 

0 - 

D = 2 C 

i - 1 

DO 

- 2€ 2 


1 

0 

» 2€ 2 

2 € 2 - 


Q 2 DQ 2 = 2C 


Using the partition (3. 27) we see that 


“1 

x 

/2 

v/Z 

X 

x 

/2 

“v/2 


= [-0 2 ] 

Qj = C/2 */2] 

Q 4 =[l] 


(3.42) 

(3.43) 

(3.44) 

(3.45) 

(3.46) 
(3. 47) 


Thus, using (3, 30) 


£l (t)= 2| 


-C(t-T) 


Vi-c 2 

a “ C (t - T ) 


-^sinVl-C' (t - t ) ' 


-3 C(t-T) 


V3-9? 


sin V 3 - 9 (t-T) 


7m? 


s in V 1 - C (t-T) 


a -3 C(t- t) 1 — 2 

;inV3- r 


-s/3 - 9 C 2 


sinV3-9C (t-T) 


(Pl (t)+P 2 (t)) 


(P 1 (7)-P 2 (T)) 


) dt 


+ 0(€*) (3.48) 


p t -C(t-T) / 2 " 

(t) = x.(t)= •—==- sinVl-C (t-T)( Pl (T)+p 2 (T)) 

3 J 0 Vl-C 2 


-2 


,-3C(t-T) 


+ sin ^9C" (t-T)( Pl (T)-p 2 (T)) 


■9C 


o -2 £(t- T) 1 o" 5 

■ r = t v - t . - sinV2-4C (t-T) Pl (T) dt+ 0(€' 

V2-4C 2 


(3.49) 
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As pointed out previously in connection with (3. 30), the response of the 
secondary system does not blow up as € tends to zero as might be suggested 
by the local rrode, whose eigenvector is [-2€, 0, € 


Example A2 [Repeated Eigenvalues] 

As a second example consider Equation (3, 0) with the matrices M, 
D and K given by 





~2 

- 

1 

0 - 


D = 2 C 

-1 


2 

-€ 2 




- 0 

- 

€ 2 

€ 2 - 



r 2 -i o- 


“0 

0 

0' 

K = 

-1 2 

0 

+ 

0 

0 

-e 2 


- 0 0 0- 


-0 

-€ 2 

€ 2 J 


(3, 50) 
(3. 51) 


(3,52) 


- 1/2 

Using the transformation xs M reduces (3. 6) to the canonical form 


PII + 

” 2 
-1 

-1 

2 

0" 

0 

+ € 

~0 0 
0 

or’ 

i 

- X 

11 

o 

(3. 53) 

X 

- 0 

0 

1. 


-0 -1 

0- 




Now, the matrix 

" 2 -l' 

-1 2 

is diagonalized by the orthogonal matrix 

~_ 1 _ jf~ 

_i_ j_ 

Jl "V2„ 

and the matrix [l] is diagonalized by the matrix [l]. Thus 
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1 

i 

0 

Vi 

Vi 

1 

1 

o 

Vi 

“Vi 



0 

i_ 


The transformation f^Tcp 1 reduces (3. 53) to the standard form 





0 

0 

i ~T 




"Vi 

IUi+ 

. 3 1. 

+ € 

0 

0 

i 

Vi 




1 

i 

o 

— 



LVi 

Vi 



(3. 54) 


(3. 55) 


We note that in this case the first and third elements of Aq are equal, this 
is then a case of repeated roots and we must use the techniques of Section 
2 applicable to the problem of repeated roots. 

First we form the auxiliary eigenequation 


[Xl + B r ]t =0 

— n 


(3. 56) 


[Note: In the formulation of Section 2 the first r_ roots were repeated; here 
the first and third roots are equal. This poses no difficulty since it can be 
put into standard form by interchanging the second and third columns of T in 
(3. 54). This is really unnecessary provided we use the elements of cV 
associated with the repeated roots. ] 

In the case of the present problem 




(3. 57) 


Solving the eigenvalue problem (3. 56), the associated orthogonal matrix T r is 
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(3. 58) 




r 


ft 

’vi 

T r = 

i 

JL 


ft 

VI_ 

The required orthogonal matrix T v is 

given 



1 

VI 

i 

VI 

1 o 

1 

VI 

0 

1 

VI 

T* = 

1 

VI 

1 

"VI 

0 

0 

l 

0 


0 

0 

1 

i 

Lvi 

0 

rS]|- w 
1 


(3. 59) 


[Note: The reordering of T r to accommodate the fact that the first and third 
roots are repeated.] Thus 




2 

1 

2 

ft 


ft 

ft 


ft 


(3,60) 


i * i 

Using the transformation jr ~ T 5,< C£ 


reduces (3. 53) to the canonical form for 


repeated roots 


W.I + 


1 0 O' 
0 3 0 
0 0 1 - 


+ € 


VI 

1 

2 

0 



0 

1 

2 


cp 1 = 0 



(3.61) 


Application of the perturbation technique applicable to repeated roots gives 


Z = 


1 

€ 

4 

0 


€ 

4 

1 

€ 

4 


0 

6 

4 

1 


(3.62) 


r 
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T 2 S t*z 

Q z = M' 1/z T*Z 


J. V* 

2 4 'Jz 

4. 

4i 

w | v» 

6 4v2 

1, € 

X 

1 6 

2 

"Vz 

2 " 4v^r 


i 


L ^ 

2\'Z 

K 


(3. 63) 


(3.64) 


The congruence matrix Q 2 has the following properties 
i) Q^MQ 2 s X + 0(6 2 ) 


T 


ii) q*kq 2 - 


l-~ 0 
V2 

0 3 


0 


0 1 + 




+ 0(6* 


(3.65) 


(3. 66) 


In this particular example we see that the first and third modes 
of vibration are both local modes and are associated with the reduction of 
the degenerate eigenvalue 1. It should be noted that in this case the relative 
motion of the secondary system is of order 1/6 compared to that of the 
main system, whereas in the previous case the relative motion was of 
order 1/6 . The second mode of vibration is a body, or global, mode which 
has eigenvalue 3 and eigenvector (^=, and is simply an extension 

of the second mode of the main structure into the secondary structure. 

Using the values of the matrix D in (3. 42), then 

€ 


,T 


Q 2 DQ 2 = 2C 


1- 


•li. 


i+ 




(3. 67) 
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Using Equation (3. 30) the solution of (3, 0) is 


*l(t) = J |*V 1 (t«T)Cp 1 (T) + P2(T)] + |-V 2 (t-T)[p 1 (T)-p 2 (T)] 

+ ~v 3 (t-T)[p 1 (r)+ P 2 (T)]dT+ 0(£) (3. 68) 

t 

x 2 (t) = f i-V 1 (t-T)Cp 1 (T)+ p 2 CT)3 - 5 -V 2 (t-T)[ Pl (T)-p 2 (T)] 

+ |-v 3 (t-T)[p 1 (T)+ p 2 (t)] dT+ 0(€) (3. 69) 


*o(t) « f [v 1 (t-T)-V 3 (t-T)]Cp 1 (T)+P 2 (T)] 

J 0 2/2 

+ ^V 2 (t-T)[ Pl (T)-p 2 (T)] 

- •|[v 1 (t-T) + V 3 {t~T)][ Pl (T)-p 2 (T)]dT+0(€) (3. 70) 


where 


V x (t) = e 


Z 1 -^) 


sm 



v 2 (t) = e 


~3£t 



(3.71) 



It will be observed that the term (t) contains a term of order l/€, One 
might think that this term will become singular as €-*0; this is not the case 
and it is easily shown that the limit of x 3 (t) as €-*0 is finite. 
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Limit (t) 
€-0 5 

— f 
" /2 0 


i_(1.2C 2 ) 

(t-T)o“' (t " T) ( ^ — * cos 

i-r 


A~C (t-T) — — ^=sin/l-c (t-T) 


~i~(l- 2 C 2 ) 

— (t-T) (p 1 (T)+p a (T)) dT 


(■-?) 

+ J |Cv 2 (t-T)-V 1 (t-T)][p l (T)-p 2 (T)] dT 


(3. 72) 


Thus we see that x^(t) is finite if C is greater than zero and and p^ 
are finite, 

Example B 

Consider now the problem 

Mx+ €Dx + Kx?= f_(t) 
x(0) = a , x(0) = b 


(3, 73) 


where M, D and K are symmetric, M is positive definite and D and K are 

* 1 - 1 

non-negative definite, Suppose that M D does not commute with M K, 
Equation (3. 73) does not possess classical normal modes. Let us first 
reduce {3, 73) to canonical form. Let 

£*M 1/2 x , cjJt)*7M” l/2 f(t) 


Then 


Iy + Cy; = gjt) 

yfO^M^a , £(0)=M 1/2 b 


(3, 74) 


Now let X”T£, where T is the orthogonal matrix which diagonalizes C. Thus 

Iz+ €Bz + A n z sr(t) 

_ 0 ' (3.75) 
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Equation (3, 75) may bo written as a system of first order equations 


where 


dw 

w(0) E3 d 


(3, 76) 


A * 




(3,77) 


w 


fl 


(3. 78) 


0 


0 


~ * ** 

0 

r.(t) 


j T a(t)_ 


o T m 

L# **«i» 


.1/2 -1 
T M. a 


Q T Ma 

T T M 1/2 b_ 


^Q T Mb^ 


(3.79) 


(3. 80) 


~l/3 

where QcM T is the congruence transformation which simultaneous iy 
diagonalizes M and K, 

The problem has now been reduced to the form where the matrix 
perturbation theory of Section 2B can be applied. In Section 2B it was 
shown that the similarity matrix T which will diagonalize A is given by 


i+j«sn 0 

x-jssn 0 ’ 

jO 0 - 6 | 

< A i +sn o) 

I j 

-jno-si 

(Aj+s^) 


(3.81) 



A x +SQq 


J - j 6 Q q 

A l +Sn 0 

jOo 1 - 53 . 


(3. 82) 
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.warn. m 


*> t 


where 


n -a 1/2 
0 " 'o 


W. 


W 


2 . 


w 


N 


(3. 83) 


where w. P u-1,2 . . . N are the radial frequencies of the undamped system. 

a. 


. _ i 
A 1 " 2 


11 




a 


22 . 


12 


2 2 
w r w 2 


®21 


“2 2 
•w 

a 1 


w 2‘ w l 


'N1 


2 2 
W N - W 1 


* a 


NN 

a 


IN 


2 2 
w -w 
1 N 


a. 


2N 


2 2 
w 0 -w_ T 
2 N 


(3.84) 


(3.85) 


Making the transformation 


Equation (3. 76) becomes 


W = t £ 


d £ 

"jn 0 -6 Aj 

0 

dt = 

0 

-3V eA J 


£+ T'bft) 




£(0) = T" 1 d 


Thus 


(3. 86) 


(3. 87) 


£(t) = 


(j° 0 - € A 1 )t 


0 e 


(-j^ 0 -^ r /C 


!J 


T ^w(0) + 


'(jQ 0 -€A 1 )(t-T) 


(“j Q o“ €A i)( fc “ T ) 


T" 1 s(T)dT (3.88) 
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Therefore 


w(t) = T 


(j V €A i ^ 


O e 




T“ 1 W(0)+ I T 

'o 


m 0 ' eA i^ 


0 o 


(-jo^eh^t 




T -1 

• 

' 0 ’ 
act) 


dT (3. 89) 


Substituting (3.78), (3,81), (3.82), (3.83), (3.84) into (3,89) and 
taking the special case where w(0)£_0 

-tl 


z (t) = 


where 


\ £ 

- ^ (8. . (t-T) sinw.(t-T) 

2 li ' ' i ' 


Q T f(T)dT 


1 \l 

t 

+ ej [sR(t-T) -R(t-T)s]Q T £ t (T)dT + 0 (€ 2 ) 


(3. 90) 


Now 


R.(T] 


\ 4 8,.t 
\ 2 11 . 
e cosw.t 


\ 


x(t) = M' 1/2 Tz(t) = Qz(t) 


l(t) = Q 
*0 


\ - 7 6,. t-T smw. t-T 
2 n ' l 

e — 

i \ 


w 


Q T l(r)dr 


A 

+ €J Q[sR(t-T)-R(t-T)s]a T _f(T)dT + 0(6* 


(3.91) 


(3. 92) 


The first term in (3. 92) is the solution which would have been obtained 
T 

if the matrix B = Q BQ had been diagonal. Thus the error in neglecting the 
off-diagonal terms is 

x(t).-x_(t) = S f Q[SR(t-T)-R(t-T)S]Q T i(T)dT + 0(€ 2 ) (3.93) 
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To obtain some idea of the errors introduced by the off-diagonal 
terms in the matrix B, let us compare the two matrices vft) and 
€[SR(t)-R(t)S], We have 



sinw jt 

e 

w l 

• 

« 



v(t) «s 

£ 

(3. 94) 



.iW sinv W 




W N __ 



and 

€[SR(t)-R(t)S] 

€ %2 (R r R 2 } € ®13 (R r R 3 } €a LN (R r R N } 

0 22 22 22 

w r w 2 w r w 2 w r w N 

€55 l2 (R r R 2 ) A €B 23 (R 2 -R 3 ) 

2 — 2 — 0 2 — 2 

w r w 2 W 2" W 3 

^l^rV n 

2 2 * U 

w r w N 


(3. 95) 


The contribution of v(t) to the solution depends on 

i) the diagonal terms 

ii) the frequencies 

iii) the frequencies and amplitudes of the model excitation 
The contribution of €[SR(t)-R(t)S ] to the solution depends on 

i) .the off-diagonal term Mj 

2 2 

ii) the frequency separation of the modes w^ -w^ 

iii) the frequencies and amplitudes of the modal excitation 
From the above one may conclude the following: 
i) The effect of the off-diagonal terms will be smaller the smaller the 
magnitude of the off-diagonal terms. 
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ii) The effect of the off-diagonal terms will be smaller the wider the 
separation in the modal frequencies. 

To illustrate some of the points in this case, consider the numerical 
Example B. 


Example B 

Consider the steady state solution of the problem 

T .., oi., r i.05 -0,051 
Iy+€ Lo 6 L-0.05 1.05 = 

The transformation 


/ cos t 
Vcos t 


JL J_ 
ft ft 

ft 


1 


AT 
v 2 


will diagonalize the stiffness matrix, giving 

Mi 1VG 


ft cos t 
0 


) 


The associated eigenvalue problem is 

(W.i IMS 

Applying the methods of Section 2B, we have 

"o - G /ft] 

*, ■ r; ;i 



o 

d~ 

i 


"o 

lo" 

II 

Lq. 1 ! 0 J 

II 

10 

0 

1 

! 1 

u J 


Q = T = 


ft ft 
_i_ J_ 

ft ’ ft 


(3.96) 


(3. 97) 

(3. 98) 

(3. 99) 
(3. 100) 

(3. 101) 
( 3 . 102 ) 
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'i’hus the steady state solution is given by (3. 92) with the lower limit set 
to -oo, 


jC(t) =x(t) =1 


A 

"i r 


f 

/z Sz 



_±_ 


J 

-00 

fZ Sz__ 



-6(t-T) . . 

e Bin (t-T) 


0 e 


-6(t~T) sin/lTl (t-T) 

/ITT 


fZ cost 


0 


df 


+ € 


j 

-00 


fz /Z 
_L J_ 
/2 "/2 


0 10[e”^" T ^][cos(t-T)-cos /IT T (t-T)] 

10[e ^ T) ][cos(t-T)-cos/n(t-T)] 0 


Jz cos t 
0 


dT + 0(6 ) 
(3. 103) 


Therefore 


X(t) = 



“* sin t~ 


- 


26 

+ 

5 cos t 


sin t 

-5 cos t 


*- 2 6 J 



+ 0 ( 6 ) 


(3. 104) 


The first term represents the response y^(t) of the system neglecting the 
off-diagonal term, thus the error due to the off-diagonal terms in B is 
given by 

^ = i(t)-Xo(t) = [Jcost] + °< 6 > < 3 - 105) 


A measure of the relative error is 

* ° sf/lfeb - 106 < 3 - 106) 

Thus we see that as long as € is small the relative error in neglecting the 
off-diagonal terms in the B matrix is small in this case. 
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